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ABSTRACT 

Spatial and spectral profiles of O VI emission behind a shock wave on the northern 
edge of the Cygnus Loop were obtained with the FUSE satellite. The velocity width 
of the narrowest O VI profile places a tight constraint on the electron-ion and ion-ion 
thermal equilibration in this 350 km s~^ collisionless shock. Unlike faster shocks in 
SN1006 and in the heliosphere, this shock brings oxygen ions and protons to within a 
factor of 2.5 of the same temperature. Comparison with other shocks suggests that shock 
speed, rather than Alfven Mach number, may control the degree of thermal equilibration. 

We combine the O VI observations with a low resolution far UV spectrum from 
HUT, an Ha image and ROSAT PSPC X-ray data to constrain the pre-shock density 
and the structure along the line-of-sight. As part of this effort, we model the effects 
of resonance scattering of O VI photons within the shocked gas and compute time- 
dependent ionization models of the X-ray emissivity. Resonance scattering affects the 
O VI intensities at the factor of 2 level, and the soft spectrum of the X-ray rim can be 
mostly attributed to departures from ionization equilibrium. The pre-shock density is 
about twice the canonical value for the Cygnus Loop X-ray emitting shocks. 
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1. Introduction 

The Cygims Loop is a nearby, bright supernova 
remnant (SNR). Beeause it is relatively unreddened, 
it has been extensively studied in the UV with lUE, 
Voyager, HUT, HST and FUSE (e.g. Raymond et 
al. 1980; Blair et al. 1991; Long et al. 1992; Van- 
cura et al. 1993; Danforth et al. 2000; Sankrit et 
al. 2000; Blair et al. 2002; Sankrit & Blair 2002). It 
is often taken as the prototypical middle-aged SNR, 
though it was probably formed by a stellar explosion 
in a relatively low density region of the interstellar 
medium (ISM) where the progenitor formed a cavity 
and a dense shell (e.g. ShuU et al. 1985; Hester et 
al. 1994; Levenson et al. 1998; Miyata & Tsunemi 
1999). In regions where the blast wave encounters 
the dense shell, radiative shocks produce bright opti- 
cal and UV emission. Regions where the blast wave 
encounters low density gas are seen in X-rays and in 
the Balmer lines as non-radiative shocks. The dis- 
tance to the Cygnus Loop is now believed to be 440 
pc based on comparison of proper motions and shock 
speeds (Shull & Hippelein 1991; Blair et al. 1999), as 
opposed to the canonical estimate of 770 pc. 

A non-radiative shock is one that heats plasma to 
a fairly high temperature, and that has encountered 
the plasma so recently that the shocked gas has had 
insufficient time to radiatively cool (McKee & Hol- 
lenbach 1980). The optical and UV emission arises 
from a narrow zone just behind the shock where the 
elements pass rapidly through successive ionization 
stages. CoUisional excitation causes some of these 
atoms and ions to emit photons before they are ion- 
ized. At visible wavelengths, non-radiative shocks ap- 
pear as faint filaments of essentially pure Balmer line 
emission (Chevalier & Raymond 1978; Chevalier, Kir- 
shner & Raymond 1980). 

Non-radiative shocks present unique diagnostics 
for shock speed and for physical processes in shocks, 
because the emission is produced before Coulomb col- 
lisions and radiative cooling erase the signatures of 
conditions at the shock. Thermal equilibration among 
different particle species is especially important, as it 
must be understood for reliable interpretation of elec- 
tron or ion temperatures in terms of the shock speed 
(Raymond 2001). 

Several non-radiative shocks along the periphery of 
the Cygnus Loop have been studied in an effort to ob- 



tain shock parameters and to investigate the physics 
of collisionless shocks. The filaments observed were 
on the western edge (Raymond et al. 1980; Treffers 
1981) and the northeastern boimdary (Raymond et al. 
1983; Fesen & Itoh 1985; Blair et al. 1999; Sankrit et 
al. 2000; Sankrit & Blair 2002). However, relatively 
low shock speeds and incipient radiative cooling made 
the interpretation of these shocks ambiguous (Long et 
al. 1992; Hester et al. 1994). Recently, Ghavamian et 
al. (2001) studied a filament on the northern bound- 
ary. It was chosen because there was no indication 
of the onset of radiative cooling. Indeed, the cooling 
time for the shock parameters derived below exceeds 
10^ years, while the flow time through the region of 
interest is 500 years. In addition, its morphology sug- 
gested a higher shock speed, making it easier to re- 
solve the Ha profile. 

Recent studies of SNR shocks have shown elec- 
tron temperatures far below proton temperatures, 
Te << Tp, in shocks faster than 1000 km s^^ (Ray- 
mond, Blair & Long 1995; Laming et al. 1996; 
Ghavamian et al. 2001) and ion temperatures ap- 
proximately mass proportional to the proton temper- 
atures, Ti ~ (mi/mp)Tp (Raymond, Blair & Long 
1995). This suggests that the plasma bulk velocity is 
randomized without efficient sharing of thermal en- 
ergy among different particle species. Slower shocks, 
in particular the Cygnus Loop shock observed here, 
are much closer to electron-ion thermal equilibration. 
Ghavamian et al. (2001) derived a shock speed Vs of 
270 to 350 km s"^ and Te/T; = 0.70-1.0 from the Ha 
line profile, the Ha to H/3 line ratio and the intensity 
ratio of broad to narrow components. Temperatures 
derived from X-ray observations in this region (De- 
courchelle et al. 1997) favor the upper ends of the 
allowed ranges, 350 km s~^ and Te/Tj = 1.0. One 
might expect the kinetic temperatures of oxygen and 
hydrogen to be close to equilibrium as well. 

A second aspect of the nature of SNR filaments is 
the 3-dimensional structure of the SNR shock. This 
is important for understanding the nature of the in- 
teraction between SNR blastwaves and density inho- 
mogencities. It is also needed to disentangle the pre- 
shock density, uq, from the depth along the line of 
sight, L, for estimates of global SNR parameters. Hes- 
ter (1987) demonstrated that SNR filaments are tan- 
gencies between the line-of-sight and a rippled sheet 
of optically thin emitting gas. The implied relation- 
ship between surface brightness and Doppler velocity 
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for the diffuse omission between the tangencies can 
be exploited to derive a 3D model of the shock (e.g., 
Raymond et al. 1988; Danforth et al. 2000). An 
alternative approach is to compare observablcs that 
depend on the line-of-sight depth, L, and the pre- 
shock density, no, in different ways. Szentgyorgyi 
et al. (2000) separated no and L in the XA region 
of the Cygnus Loop by comparing the intensity of a 
[Ne V] line, which scales as noL, with the X-ray sur- 
face brightness, which scales as tiqL. Patnaudc ct 
al. (2002) applied a similar technique to a recently 
shocked cloud in the SW Cygnus Loop. 

A quantitative study of the O VI emission requires 
consideration of the effects of resonance scattering, 
both within the emitting sheet of gas and in the in- 
tervening ISM. The intrinsic intensity ratio of the res- 
onance doublet of a Li-likc ion is 2:1, but the 2:1 ra- 
tio of opacities implies a larger optical depth for the 
shorter wavelength line. Thus departures from a 2:1 
ratio signal non-negligible optical depths in the O VI 
lines (Long et al. 1992; Sankrit & Blair 2002). Cor- 
nett et al. (1992) compared the optically thin emis- 
sion in [O III] A5007 with UV omission dominated 
by the C IV A1550 doublet to show that bright fila- 
ments in the eastern Cygnus Loop are optically thick 
in the C IV lines. They confirmed that the filaments 
are tangencies to the line of sight and that resonance 
scattering is important for strong UV lines. 

This paper presents O VI line profiles obtained 
with the FUSE satellite at 4 positions behind the fil- 
ament investigated by Ghavamian et al. (2001) along 
with a low resolution FUV spectrum from the Hop- 
kins Ultraviolet Telescope (HUT). The line widths 
are used to determine the degree of ion-ion thermal 
equilibration. The spectra yield both the line pro- 
file information and the intensity ratio information 
to permit both of the approaches described above to 
constraining the 3D structure of the shock. We con- 
struct models of the emission and scattering in the 
sheet of shocked gas to interpret the spectra. We find 
nearly complete equilibration between O and H, a pre- 
shock density about twice the value usually quoted 
for the Cygnus Loop blastwave, and a depth along 
the line of sight of about 0.7-1.5 pc. We also ana- 
lyze the ROSAT PSPC X-ray imago of the region. A 
thin, soft rim appears behind the filament in ROSAT 
data. We construct time-dependent ionization models 
to resolve the ambiguity between lower shock speed 
and emission from lower ionization states as causes of 



the soft rim. Wo find that the soft rim at this loca- 
tion is largely due to strong emission from moderate 
ionization states in the region close to the shock. 

2. Observations and Data Reduction 
2.1. FUSE Observations 

Spectra were obtained on 15 and 16 June, 2000 
with the Far Ultraviolet Spectroscopic Explorer satel- 
lite (FUSE; Moos et al. 2000, Sahnow et al. 2000). A 
detailed discussion of the procedures and uncertain- 
ties related to FUSE observations of SNR filaments 
is given by Blair et al. (2002). Pointings were made 
at the four positions listed in Table 1 and shown in 
Figure 1. The uncertainty in the FUSE blind off- 
set is believed to be about 1". The image in Figure 
1 is a superposition of Ha (red) and [O III] A5007 
(green) images obtained through narrowband filters 
on the 1.2 m telescope at the Fred Lawrence Whipple 
Observatory on Mt. Hopkins on Nov. 7, 2000. Condi- 
tions were not photometric, and the seeing was about 
1.6". The image shows a portion of the Ha filament 
complex that defines the northern edge of the Cygnus 
Loop. Faint Ha emission well out in front of the main 
filaments reveals a portion of the blastwave in lower 
density gas. This is also apparent in the ROSAT im- 
age as faint X-ray emission ahead of the main edge 
of the remnant. A portion of the shock making the 
transition from non-radiative to radiative appears in 
[O HI] emission at the southeast corner of the Figure 
1. Very faint, diffuse [O II] and [O III] emission in 
the northeastern part of the image may arise from the 
partially ionized pre-shock gas (Bohigas et al. 1999). 

As shown in the inset in Figure 1, the 20" by 4 " 
MDRS aperture was placed parallel to the Ha fila- 
ment, at a position angle of 317° E of N. The aper- 
ture was placed on the brightest Ha filament at the 
same position observed in the optical by Ghavamian 
et al. (2001). Subsequent slit positions were parallel 
to the first, but displaced 6.5", 11", and 15.5" be- 
hind the first. The last position falls upon a fainter 
Ha filament located about 15" behind the bright one. 
Because of imccrtain alignment between the LiF and 
SiC channels, we consider only the LiF data for the 
MDRS spectra. 

Spectra through the 30" by 30" LWRS aperture 
were obtained simultaneously. The LWRS aperture 
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was located about 3.5' to the northwest along the 
same filament complex, at the positions shown in Fig- 
ure 1. The LWRS sensitivity is higher by the ratio of 
aperture areas. Moreover, the SiC channels can be 
used to extend the spectral coverage because the oflf- 
set between LiF and SiC channels is relatively small 
compared to the LWRS aperture size. However, the 
LWRS significantly degrades the spectral resolution 
for extended emission. Especially when optical depth 
within the emitting gas or in the intervening ISM af- 
fects the line intensities, it is more difiicult to de- 
rive reliable estimates of the intrinsic emission inten- 
sities. Therefore, the LWRS spectra arc less useful 
for the purposes of this paper, and we discuss them 
only briefly below. 

The data were processed with CalFUSE 2.0. We 
extracted MDRS spectra from the complete data set 
and also extracted spectra from the night portions of 
the orbits alone. The night-only extraction greatly 
reduced airglow, but it reduced the signal-to-noise of 
the O VI profiles so severely that we present only 
the data from the full exposures. For each MDRS 
pointing we shifted the Lifla and Lif2B spectra to 
align the O I A1039.4 airglow features, then added 
the Lifla and Lif2b spectra weighted by the effective 
areas of those channels. The residuals around the 
fit to the dispersion solution for the LiF 2B segment 
are around 5 km s~^. The Lif2B aperture may have 
drifted relative to the LiflA aperture by a fraction of 
the slit width during the course of the observation, 
effectively broadening the aperture by 1 or 2" (Blair 
et al. 2002). The absolute radiometric calibration of 
FUSE is believed to be good to 10% , and the relative 
calibration between 1032 A and 1038 A should be even 
better. 

Figures 2a and 2b show the profiles of the O VI 
1032.93 A and 1037.61 A lines at the four positions, 
with the uppermost profiles corresponding to the posi- 
tions closest to the blastwavc. Interstellar absorption 
in the C II A1037.02 and H2 AA1037.15, 1038.16 lines 
strongly affects the 1037.61 A line. While the inter- 
stellar lines lie at velocities of -171, -134 and -t-158 
km s""'^ with respect to the O VI line, we have no 
means to assess their effects except to compare the 
1032.93 A and 1037.61 A fine profiles. We assume 
that these absorption lines are important wherever 
the ratio I(1037)/I(1032) falls below its intrinsic value 
of 0.5. On that basis, interstellar absorption signifi- 
cantly attenuates the 1037.61 A emission at velocities 



below about -60 km s and above about -|-90 km s . 

Figure 3 shows the profile of O VI A1032 from 
the sum of the LWRS exposures farther up the fil- 
ament. The lower resolution of the LWRS smears the 
profile enough to mask the double peak structure if 
it is present. The asymmetry suggests that the dif- 
fuse region within the LWRS apertures is redshifted. 
The profile is similar to what one would expect from 
adding the profiles in Figure 2a and degrading the 
resolution to that of the LWRS, though the situation 
is more complicated because the intensity distribu- 
tion across the LWRS aperture also affects the pro- 
file. One novel feature is present, however. The faint 
wing at redshifts above 200 km s^^ probably repre- 
sents an additional component along the line of sight. 
Presumably this is due to the shock front associated 
with the faint Ha and X-ray emission ahead of the fil- 
ament observed by FUSE. If it's speed is 350 km s~^ 
like that of the filament observed by FUSE, the shock 
must be viewed at about 35° to give the observed 
Doppler shift. 

Table 2 presents the fluxes of the O VI lines at the 
four positions. It should be kept in mind that the 1037 
line is attenuated by H2 and C II interstellar lines, 
and that the effect of this absorption depends on the 
intrinsic line profile. The dereddened fluxes are based 
on E(B-V) = 0.08 (Miller 1972) and the extinction 
curve of Cardelli, Clayton & Mathis (1989). Ha fluxes 
computed by scaling the Ha brightness in each MDRS 
aperture in the image in Figure 1 to surface brightness 
of a long slit spectrum of Ghavamian et al. (2001) 
are given for comparison. This scaling should lead to 
fluxes accurate to about 20%. 

For completeness. Table 2 includes the O VI and 
C III A977 fluxes measured from the sum of the 4 
LWRS spectra. Because of the higher sensitivity of 
the LWRS spectra and the summation of 4 spectra, 
we avoid airglow contamination by using only data 
obtained during orbital night. The C HI A977 line is 
detected, but variations in the background subtrac- 
tion limit the accuracy of the measured flux to ±40% . 
Upper limits for other lines, such as Nc V] 1136.5 and 
Ne VI] 1005.7 are 2% the intensity of O VI A1032. Up- 
per limits on the S VI AA933.4, 944.5 lines are about 
6% of O VI A 1032 because of the lower sensitivity at 
the shorter wavelengths. The C III line is badly af- 
fected by interstellar absorption, so it is diflacult to 
interpret. If resonant scattering attenuates the O VI 
intensity by a factor of 2 (see section 3.3.3), the strin- 
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gent limits on No V] and No VI] as compared to radia- 
tive shock models (Hartigan, Raymond & Hartmann 
1987) confirm that this is a non-radiative shock. 

2.2. HUT Observations 

The Hopkins Ultraviolet Telescope (HUT; David- 
sen et al. 1992) observed the non-radiative shock dur- 
ing the Astro-2 mission during March 1995. The 10" 
by 56" slit was centered at the position of the first 
FUSE observation with a position angle of 315°. The 
HUT aperture encompasses FUSE positions 1 and 2, 
and it extends nearly 3 times as far along the filament. 
Essentially all of the time when the Cygnus Loop was 
observable fell during orbital day, so the geocoronal 
background was high. 

The HUT spectra cover the 912-1800 A range with 
3-4 A resolution, and the radiometric calibration is 
accurate to about 10% (Kruk et al. 1999). Only the 
hues of He II, C IV and O VI were detected. The 2- 
photon continuum of hydrogen, C HI A977 and the N 
V doublet at 1238, 1242 A are expected to be strong 
in non-radiative shocks, but we could not reliably sep- 
arate these features from the wings of geocoronal Lya 
or Ly7. 

Figure 4 shows the wavelength range that covers 
the C IV and He II lines as extracted from expo- 
sures totahng 2861 seconds. The O VI 1032 and 1037 
lines are badly compromised by geocoronal emission 
at this resolution. Therefore, we report only the in- 
tensities of C IV AA1548, 1550 and He II A1640. Com- 
parison of measurements from the first and second 
pointings (1311 and 1550 seconds, respectively) indi- 
cates that fluctuations in the background level and 
the uncertainty in the double Gaussian fit dominate 
the uncertainties. We find intensities 3.20±0.6 and 
4.38±0.8 X 10~^^ erg cm~^ s~^, respectively. 

The lack of intercombination lines such as N IV] 
A1786 and O IV] A1400 confirms the non-radiative 
nature of the shock. In the models of radiative shocks 
these lines are comparable to C IV and He II in bright- 
ness (Hartigan, Raymond & Hartmann 1987). The 
lack of O IV] provides the stronger constraint, demon- 
strating that the shocked gas has not yet cooled to 
2 X 10^ K. 



2.3. ROSAT Observations 

ROSAT PSPC observations of the Cygnus Loop 
have been presented by Levenson et al. (1999). The 
filament observed here is about 20' east the center of 
their 3506 second North (E) region exposure. For the 
purposes of this paper, we extracted PSPC spectra 
of fourteen 25" by 250" bands parallel to the fila- 
ment for comparison with theoretical shock models. 
The 25" width is comparable to the spatial resolu- 
tion of the PSPC; so the spectra are not entirely in- 
dependent. Single temperature fits with a Raymond 
& Smith (1977) ionization equilibrium model yield 
temperatures increasing from 0.14 to 0.2 keV over a 
distance of about 1' behind the shock. However, ion- 
ization equilibrium is a poor approximation close to 
the shock. We therefore fold the spectra predicted by 
the shock models through the ROSAT effective area 
to compute the PSPC count rate as a function of po- 
sition behind the shock. The models and results are 
discussed in section 3.3.4. 

3. Analysis 

3.1. Resonance Scattering 

Before attempting to interpret the line profiles, we 
need an estimate of the optical depths in the O VI 
lines. Models based on the code described by Ray- 
mond (1979) and Cox & Raymond (1985) predict op- 
tical depths of 0.01 to 0.1 along the flow direction. 
The filament observed by FUSE is viewed nearly edge- 
on, so the optical depths are much larger. Because 
the emitting sheet of gas is much thicker along the 
line of sight than in the flow direction, the main ef- 
fect of finite optical depth is to scatter photons out 
of the line of sight. As the optical depth in the A 1032 
line is twice as large as that in A1037, the intensity 
of A1032 is more strongly affected, and the intensity 
ratio I(1032)/I(1037) drops below its intrinsic value 
of 2:1 

For a slab geometry and the approximation of sin- 
gle scattering, the intensity of a line is given by 

7 = e (1 - e--) (1) 

where e is the ratio of emissivity to opacity. The 
factor e includes the collision strength, the oscillator 
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strength, the electron density and factors involving 
the electron temperature and the O VI line width. 
All these factors are nearly constant in the ionization 
zone just behind a non-radiative shock, so e is nearly 
constant. In that case, the effects of resonance scat- 
tering can be assessed from the ratio 



R 



2 7(1037) - 7(1032) 2(1 - e-'"i°37) - (1 - e-''^" 



2 7(1037) 



2(1 - e-^ios'-) 
(2) 

with T1032 = 2tio37. The ratio varies from for r=0 
to 0.5 for T=oo. 

Figure 5 shows the observed ratios R for Positions 1 
and 4, and Figure 6 shows the function given in Equa- 
tion 2. The ratios for Positions 2 and 3 are roughly 
0.1, but they are quite noisy. The solid parts of the 
curves are the velocity ranges where the data are re- 
liable. The dotted portions are regions where C II or 
H2 interstellar absorption affects the profile, or where 
the signal in the 1037 line is too low to give a reliable 
ratio. 

Figures 5 and 6 imply optical depths rio37 of or- 
der 1-3 through the emitting regions of Positions 1 

and 4. At Positions 2 and 3 the ratios are smaller, 
suggesting optical depths of 0.3 to 0.5. The smaller 
optical depths are consistent with the morphological 
indication that the line of sight is farther from tan- 
gency at positions 2 and 3. In a small portion of 
the range in Figure 5a near a velocity of km s~^ R 
exceeds the theoretical maximum of 0.5. This is pre- 
sumably due to interstellar absorption, which scatters 
A1032 photons more effectively than A1037 photons, 
but does not compensate by emitting correspondingly 
more A1032 photons. However, Figure 5b does not 
show a corresponding peak at position 4. This sug- 
gests that a combination of intrinsic and interstellar 
scattering is required to account for the peak in Fig- 
ure 5a. 



3.2. Thermal Equilibration 

The discussion above shows that optical depth can 
affect the line profiles, but the rough constancy of R 
in Figure 5 suggests that scattering affects the overall 
intensity more than it distorts the line profile. Com- 
parison of the O VI profiles obtained from the the 



LiF lA and LiF 2B c;liannc;ls indicates that the fluc- 
tuations in Figure 5 are marginally real. The most 
stringent limit on the line width comes from the sharp 
feature a.tV= +25 km s~^ in the profile for Position 
4 that arises from the trailing filament seen in Figure 
1. Based on Figures 5b and 6, the optical depth at 
the peak of the bright feature in the A1032 profile at 
Position 4 is about 1.5, and the flux is attenuated by 
about a factor (1 — e~'^)/T, or about 0.52. The fluxes 
at the apparent half-power points at -20 km s"^ and 
+60 km s~^, on the other hand, are attenuated by fac- 
tors of 0.63 and 1.0, respectively. Therefore, the line 
profile has been broadened by attenuation at its peak, 
and the apparent Full Width Half Maximum of 120 
km s~^ corresponds to a 90 km s~^ intrinsic width. 
It is clear from the profiles of Positions 1 to 3 and 
from the blue wing of the Position 4 profile that bulk 
motions also contribute to the Doppler width. There 
is no reliable way to assess the bulk velocity contri- 
bution, however, so we use the observed line width 
corrected for resonance scattering as an upper limit 
to the oxygen kinetic temperature 



To < 2.7 X 10*^ K 



(3) 



This can be compared with the proton temperature 
of Tp = 1.4 ± 0.35 X 10^ K derived from the Ha line 
width by Ghavamian et al. (2001). The time scale for 
0^+ ions to slow by Coulomb collisions corresponds 
to a distance of ~ 10-15" behind the shock for the 
density range derived below (Spitzer 1968), while the 
O VI emissivity peaks only 1" behind the shock in the 
model presented below. Therefore, the upper limit in 
To applies to the immediate post-shock region. Thus 
the upper limit to the ratio of oxygen to hydrogen 
kinetic temperatures is 2.5 including the uncertainty 
in the proton temperature. 

For comparison, Raymond et al. (2000) and Man- 
cuso et al. (2002) found O VI line widths near 3/4 
Ys (FWHM) for 1000 km s"! shocks in the solar 
corona. The 350 km s~^ Cygnus Loop shock shows 
a line width below 1/4 Vg. The Alfven Mach num- 
ber, defined as the shock speed divided by the Alfven 
speed, might be a parameter that controls heating in 
the shock. In spite of their high speeds, the shocks 
in the solar corona had modest Alfven Mach num- 
bers and modest compressions. The Cygnus Loop 
shock has a lower speed, but probably a higher Alfven 
Mach number. Berdichevsky et al. (1997) found 
He and O kinetic temperatures roughly twice the 
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mass ratios times the proton temperatures in shoeks 
in the solar wind. Few observations of ion temper- 
atures in supernova remnant shocks are available. 
Roughly mass-proportional temperatures, Tq ^ I'iTp 
and Tffe ~ 4Tp, were found in a 2300 km s~^ shock 
in SN1006 (Raymond, Blair & Long 1995). 

We find that for this 350 km s'^ shock To Tp - 
Te even though its Alfven Mach number is fairly high, 
as opposed to faster shocks in which Tq > Tp > Tg. 
On the other hand, fast shocks in the solar wind with 
modest Alfven Mach numbers also show To > Tp > 
Te. This suggests that shock speed rather than Alfven 
Mach number may control the thermal equilibration. 

3.3. Structure in 3 Dimensions 

It is not practical to invert the observed Ha image 
and the O VI profiles to obtain the shape of the shock, 
because the effects of resonance scattering are compli- 
cated and sensitive to both the viewing angle and the 
geometrical structure (e.g., Wood & Raymond 2000). 
Therefore, we proceed by constructing a model, com- 
puting the line intensities and profiles for the 4 po- 
sitions, and adjusting the model parameters to opti- 
mize the fit. We first discuss the overall nature of the 
structure and its relationship to the measured pro- 
files, then describe the model, and finally present the 
results. 

3.3.1. Model Framework 

The Ha filaments clearly represent the projection 
of a rippled sheet on to the plane of the sky as pro- 
posed by Hester (1987). The surface brightness of the 
shock seen face-on is the product of the number of H 
atoms swept up per second and the number of pho- 
tons emitted by an H atom before it is ionized, which 
is given by the ratio of excitation and ionization rates. 

Iq = —— TiHo photons cm~^ s~^ sr~^ (4) 

47r qion 

For the nearly edge-on viewing angle of the filaments, 
the brightness is increased by a factor sec(6'), where 
9 is the angle between the line-of-sight and the shock 
normal. As 6 approaches 90°, the geometrical en- 
hancement is limited to L/lobs, where L is the depth 



along the line-of-sight and lobs is the larger of the 
thickness of the emitting layer and the spatial resolu- 
tion element. The two nearly parallel bright filaments 
covered by the FUSE MDRS slit positions are tangen- 
cies between the line-of-sight and the emitting sheet, 
and the diffuse region in between is where the sheet 
curves away from tangency. 

An equation similar to Equation 4 describes the 
O VI intrinsic brightness, but with no replacing uho 
and the excitation and ionization rates for the O VI 
lines replacing those of Ha. However, the projection 
of the O VI brightness is more complicated. The 
thickness of the Ha emitting sheet is smaller than the 
resolution of the images, < 10^^ cm in the models 
presented here, and the Ha line is optically thin. For 
the O VI lines, on the other hand, the emitting sheet is 
roughly 5" thick. The optical depths discussed above 
present an even greater complication. The line pro- 
files are determined by the emissivity as a function of 
position, the line-of-sight component of the bulk ve- 
locity of the shocked gas, the thermal velocity profile 
of the O VI ions, and the optical depth as a function of 
position and velocity. We adopt an oxygen abundance 
of 8.73 from Holweger (2001). This is lower than the 
widely used value of 8.82, but it is supported by the 
work of Allende Prieto, Lambert & Asplund (2001). 
The emitting layer is thin enough that radiative cool- 
ing has no effect, so the abundances of other elements 
do not affect the predictions for O VI. 

3.3.2. VI Profile Models 

The emissivity and opacity behind a planar shock 
wave arc computed by the code described in Raymond 
(1979) and Cox & Raymond (1985), which has been 
updated to allow electron and ion temperatures to 
evolve separately. The model parameters are shock 
speed, pre-shock density, pre-shock ionization state, 
pre-shock magnetic field, and the ratio of electron to 
ion temperature at the shock. The magnetic field 
and pre-shock ionization state have little effect on 
the O VI emission and absorption, and they are fixed 
at a neutral hydrogen fraction of 0.5 and a field of 
3 fiG. The electron-proton temperature equilibration 
was investigated by Ghavamian et al. (2001), and 
in the Te/Tp > 0.7 range they found, the initial 
electron-ion temperature ratio also has little effect. 
In the extreme non-equilibrium case Tg/Tp = 0.05, 
the thickness of the O VI emitting layer is roughly 
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twice as large as it is for the Te/Tp = 1 models we 
adopt. 

The model code starts with the Rankine-Hugoniot 
jump conditions, then computes the density, velocity, 
Te and Ti, and the time-dependent ionization state 
of the shocked gas. Because we are interested only 
in the narrow ionization zone behind the shock, the 
density is very close to 4 times the pre-shock density, 
and the velocity in the shock frame is very close to 
1/4 the shock speed throughout. For nearly complete 
electron-ion equilibration in the shock, the temper- 
ature is also nearly constant in the emitting region. 
Thus the only important parameters are Vg and ng. 
The O VI excitation and ionization rates are fairly 
weak functions of temperature above 10^ K, and the 

1 /2 

O VI thermal line width scales as TJ . The thick- 
ness of the shock structure scales as , the opacity 
scales as uq and the local emissivity as Uq. 

The shape of the shock front along the line of sight 
is parameterized as 

Y = a X -I- b sin(7rX/c) (5) 

where X is the distance along the line of sight and Y 
is distance in the plane of the sky. While this form 
imposes a symmetry on the structure, the data do 
not warrant a function having more parameters. The 
models used a pixel size of 6.6 x 10^'' cm, or 1.0" at 
440 pc. We also assume that the shock speed and 
preshock density are constant along the shock front. 
Allowing them to vary would introduce unresolvable 
ambiguities. We use the density, temperature and 
O VI concentration from the planar shock models to 
compute the opacity and emissivity in the O VI lines 
at positions behind the shock. The opacity and emis- 
sivity are computed at 10 km s^^ resolution using 
the O VI line width, shifted by the line-of-sight com- 
ponent of the bulk speed. We adopt an oscillator 
strength fio32 — 0.13 (Morton 1991), an excitation 
rate coefficient qio32 = 1-5 x 10~* cm^ s~^ (Zhang, 
Sampson & Pontes 1990), and an ionization rate from 
Younger (1981) . 

The interstellar O VI absorption is an "external" 
parameter of the models. Shelton & Cox (1994) found 
that the galactic O VI column includes a Local Bubble 
component of Nqvi = 1-6 x 10^^ cm~^ and absorbing 
regions of Nqvi = 2 — 7 x 10^^ cm~^ separated by 
400-1300 pc. The Cygnus Loop emission will be ab- 
sorbed by the Local Bubble component and perhaps 
by one other absorbing region. The ISM absorption 



is taken to have a temperature of 10^' K. The widths 
of individual O VI absorption components may corre- 
spond to 3 X 10^, but multiple components along the 
line of sight can result in typical line widths corre- 
sponding to lO'^ K (Jenkins & Meloy 1974). The Lo- 
cal Bubble alone produces a dip of about 15% at zero 
velocity. The central peak in Figure 5a suggests some 
interstellar absorption in that it rises above 0.5. It is 
consistent with a total column Nqvi = 7xl0^^cm~^. 
The model code computes the surface brightness 

I{Y,V) = ^ [ e{X,Y,V)e-'''^^'^'^Ux photons cm'^ s"^ sr"^ 
47r J 

(6) 

for the A1032 and A1037 lines, and we average the 
profiles over 4" in Y to compare with the FUSE ob- 
servations. Figure 8 shows model predictions for the 
A 1032 profiles at the 4 slit positions in the same for- 
mat as the observed profiles in Figure 2. As in Ta- 
ble 2, we assume a factor of 2.8 attenuation due to 
reddening. This model assumes Vs = 350 km s"\ a 
pre-shock density of 0.5 cm~'^ and an interstellar col- 
umn density Nqvi = 0.7 x 10^^ cm"^. The shape 
parameters of the emitting sheet were a = -0.17, b 
= 28 and c = 300. The shape is shown in Figure 7, 
but note the order of magnitude difference in the X 
and Y scales. This model assumed an oxygen kinetic 
temperature equal to the proton temperature in the 
O VI emitting region. Figure 9 shows the predictions 
for the intensity ratio R shown in Figure 5. 

To examine the effect of changes in the pre-shock 
density, we computed the O VI A 1032 profiles for a 
model with no = 0.3 cm~^. The total number of pho- 
tons emitted scales directly as no, and the profiles in 
Figure 9 are very close to those in Figure 7a scaled by 
a factor of 3/5. The optical depths in both models are 
similar, because the density in the no = 0.3 cm^'^ is 
smaller, but the thickness of the O VI layer is larger by 
the same factor because the cooling length increases. 
There are departures of order 10% from the 3/5 scal- 
ing because the areas extracted for FUSE positions 1 
through 4 correspond to slightly different parts of the 
post-shock ionization structure. 

To illustrate the effects of optical depth. Figure 10 
shows the intrinsic profile emitted at position 1 and 
the profile after the c;fFects of resonant scattering and 
interstellar absorption. The observed profile is atten- 
uated by a further factor of 2.8 by reddening. The dip 
at line center is largely due to interstellar absorption. 
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The profile without ISM absorption would bo nearly 
flat between the two peaks. The net effect of scat- 
tering, which is dominated by scattering within the 
emitting sheet, is to reduce the O VI flux by about a 
factor of 2. 

The resonance scattering will also affect the C IV 
intensity measured by HUT. However, the effect may 
be less severe than for the FUSE observations, be- 
cause the HUT aperture is larger and more of it is 
filled by diffuse emission that arises from more face- 
on portions of the shock. The relative intensities of 
the C IV and He II lines are in reasonable agreement 
with the predictions of models such as that described 
in section 3.3.1 for a 350 km s~^ shock provided that 
carbon is depleted by a factor of 2-3 and that reso- 
nance scattering (perhaps dominated by the interstel- 
lar C IV absorption) further attenuates the C IV flux 
by a factor of 2-3. 

As a check on the model parameters we compare 
the Ha fluxes at 4 FUSE positions with the figures in 
Table 2. Each neutral hydrogen atom passing through 
the shock emits 0.25 Ha photons on average, and the 
shock emits 

Io{Ha) = 0.25 Uho V;/47r photons/cm^ s sr (7) 

and the area within each FUSE position can be com- 
puted from Equation 5. The computed Ha fluxes are 
0.019, 0.017, 0.022 and 0.025 photons cm-^ s"^ for 
4" by 20 " boxes at the four positions respectively for 
Uho = 1.0. The fluxes at the flrst and last positions 
are extremely sensitive; to the exact position, chang- 
ing by nearly a factor of 2 with shifts of 1" in the po- 
sition. The variation with position roughly matches 
that in Table 2, but with significant scatter. The 
ratio of the total of the predicted fluxes to the to- 
tal of the observed Ha fluxes is 4.0 uho. Thus the 
predictions match observations for no=0.5 and a pre- 
shock neutral fraction of 0.5. Calculation of the flux 
of ionizing photons from such a shock, mainly He*^ 
A584 and He"*" A304, suggests an upper Hmit on the 
pre-shock neutral fraction of 0.7 (e.g.. Smith et al. 
1994; Ghavamian et al. 2000). Comparison of the 
O VI and Ha fluxes in Table 2 with the model pre- 
diction I1032/ Ihcx ~ 12 (and considering the factor of 
2 attenuation of /1032 derived above) suggests a neu- 
tral fraction near 1. On the other hand, Ghavamian 
(1999) used the ratio of the He II A4686 line to H/3 
in this fllament to derive an upper limit of 0.2 for 
the pre-shock neutral fraction. This would require 



a preshock density uq > 1.25. Such a high density 
would imply an extremely small depth along the line 
of sight to match the ROSAT X-ray brightness (see 
section 3.3.4), of order 0.1 pc. This small depth would 
increase the curvature of the sheet by a factor of 7, 
resulting in a very large separation of the peaks in the 
O VI profile. Thus wo prefer the parameters no ^ 0.5 
and a high neutral fraction. 

3.3.3. Comparison to Observed O VI Profiles 

Figure 8 shows overall qualitative agreement with 
the intensity and a line shape that evolves from double 
peaked near the primary shock toward single peaked 
where the emission from the trailing tangency domi- 
nates. As in the observed spectra, the flux is highest 
at position 1, declines through positions 2 and 3, then 
increases again at position 4. 

In detail, there arc two discrepancies. First, the 
separation of the peaks at position 1 is predicted to be 
only about 85 km s~^, while the observed separation 
is 110 km s"-'^. The predicted separations increase to 
95 and 105 km s~^ at positions 2 and 3, closer to 
the observed separations, but still somewhat smaller. 
Second, while the predicted red peak at position 4 is 
closer to zero velocity than at position 3, as observed, 
the predicted blue peak is nearly as strong as the red 
peak, while the observed blue peak is much fainter. 
The faint blue peak probably results from lower pre- 
shock density in the region producing that part of the 
profile. The separations of the peaks at positions 1-3 
present a more substantial problem. 

The separations of the peaks are partly due to the 
interstellar absorption and partly to the linc-of-sight 
velocity components of the shocked gas. The ISM ab- 
sorption cannot plausibly be made much broader, but 
a higher shock speed or a larger curvature, yielding 
larger angles between the line of sight and the shock 
front, would increase the separation. A higher cur- 
vature would require a compensating density increase 
to maintain the surface brightness, and no = 0.5. is 
already rather high in comparison with X-ray obser- 
vations discussed below. A higher shock speed would 
be difficult to reconcile with the Ha line width mea- 
sured by Ghavamian et al. (2001). 

The predicted separations of the peaks could also 
be increased by relaxing some of the assumptions in- 
herent to the model. We have assumed the flow to be 
perpendicular to the shock surface, as will occur for 
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a shock driven by gas pressure. An oblique magnetic 
field can introduce a velocity component along the 
shock surface, but for the parameters considered here 
and fields below 10 /iG, this will less than about 1 
km s~^ according to the shock jump conditions (Mc- 
Kee & Hollenbach 1980). It is also possible to main- 
tain oblique shocks if the ram pressure of the flow is 
significant, for instance as occurs in bow shocks. The 
steepest angle in the model sheet is 6°, so the paral- 
lel component of the post-shock flow speed should be 
small, but flows along the shock at ~ 20 km s~^ in 
the sense expected for oblique shocks could account 
for the observed peak separations. 

Another explanation would involve cosmic ray ac- 
celeration. We have assumed the the post-shock flow 
speed of a gasdynamic shock, 3/4 Vs. If a signif- 
icant fraction of the energy dissipated in the shock 
goes into accelerating cosmic rays, the shock compres- 
sion is larger than 4, and the post-shock flow speed is 
larger than 3/4 Vs. This explanation would require 
higher shock speeds to match observed Ha profile and 
X-ray temperature. If most of the energy dissipated 
goes into cosmic rays, a shock several times faster may 
be needed (Boulares & Cox 1988). The higher shock 
speeds would imply a greater distance to the Cygnus 
Loop in order to match measured filament proper mo- 
tions. The O VI and Ha surface brightnesses are pro- 
portional to the number of particles swept up per sec- 
ond, so the higher shock speed would imply a smaller 
density. 

3.34. ROSATData 

It is important to check the consistency of the de- 
rived parameters against the ROSAT X-ray images of 
the Cygnus Loop. Levenson et al. (1999) included 
the region we observe in their study of the narrow 
soft rim in PSPC images. Miyata & Tsunemi (1999) 
observed an area just to the East in a study of el- 
emental abundances and they favored idea that the 
blastwave has recently reached the wall of a cavity. 
Decourchelle et al. (1997) studied the region just to 
the West, exploring the multiple temperature compo- 
nents needed to explain the brighter, softer emission 
regions. 

Figure 11 shows the ROSAT counts extracted in 
25" by 250" bins parallel to the filament. The po- 
sition of the main Ha filament is taken as zero, and 
440 pc is taken as the distance to the Cygnus Loop 



(Blair et al. 1999). Models were computed with 
the Raymond (1979) code and Allen (1973) abun- 
dances. The curves show model predictions for a 
350 km s~^ shock with preshock densities of 0.2 and 
0.3 cm~'^. The X-ray emissivities predicted by the 
shock model are attenuated by an interstellar column 
Nh = 1.5 X 1020 cm-2 and multiplied by the ROSAT 
PSPC effective area and exposure time. A somewhat 
higher Nh is often assumed for the Cygnus Loop, but 
Decourchelle et al. (1997) require Nh ~ 1 x 10^° cm"^ 
to match the very soft X-ray spectra just to the west 
of the position we observed. The model predictions 
have been added to an assumed background Cygnus 
Loop emission of 1500 counts estimated from the 
brightness ahead of the filament. The predicted X- 
ray counts were scaled to line-of-sight depths of 2 and 
4 pc, respectively, for no = 0.3 and 0.2 cm~'^. The 
higher pre-shock density no = 0.5 used for the O VI 
profile models would produce a curve higher and nar- 
rower than the solid curve in Figure 11, but the spa- 
tial resolution of the PSPC would smear it to agree 
reasonably well with the measured values if the line 
of sight depth is scaled to about 0.7 pc. The PSPC 
resolution would smear the 0.2 cm~^ predictions to 
be broader than the observed points, so we conclude 
that no > 0.2cm~^. 

The post-shock density and temperature are essen- 
tially constant, and the peak in the predicted bright- 
ness results from the presence of moderate ionization 
stages that emit efficiently in the soft (0.25 keV) part 
of the ROSAT bandpass before the gas reaches ion- 
ization equilibrium. This departure from ionization 
equilibrium also produces a lower temperature in sim- 
ple fits to the ROSAT spectrum, accounting for the 
gradual rise from kT = 0.14 keV to kT = 0.2 keV 
over a distance of about 5 x lO^'' cm in the fitted 
temperatures. 

The models are planar, so any bending of the 
surface will broaden the profile compared with the 
model profiles. The 0.3-0.5 cm~^ models require a 
path length comparable to that suggested by the O 
VI observations, while a preshock density of 0.2 re- 
quires twice that value. A lower interstellar column, 
Nh = 10^" cm^2 would bring the no = 0.2 model 
into agreement with the O VI model and push the 
higher density model to a shorter length scale. We 
conclude that no = 0.3 — 0.5cm^'^ and L=0.7-1.5 pc 
arc consistent with both the FUSE and ROSAT data. 
We note, however, that the soft X-ray rim predicted 
by the model arises mostly from C, Si and Fe emis- 
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sion lines, and that those elements may be locked up 
in grains, as suggested by the C IV flux measured by 
HUT. If they are liberated by sputtering over a swept- 
up column of order 10^* cm~^ (Vancura et al. 1994), 
the peak will be suppressed. 

The soft rim seen in the ROSAT hardness ratio 
map of the Cygnus Loop was used by Levenson et 
al. (1999) to support the idea that the Cygnus Loop 
blastwave recently encountered the shell of a cavity 
created by the progenitor star. As Levenson et al. 
(1999) point out, the soft emission in the X-ray bright 
regions can result from a slower shock speed in denser 
material. They note that in the non-radiative regions 
such as the one examined here, the rim could result 
either from a higher density when the shock encoun- 
ters the shell, or from the non-equilibrium ionization 
effect illustrated in Figure 11. We find that non- 
equilibrium ionization can account for the soft rim 
in this region, but only if the refractory elements are 
present in the gas phase without serious depletion. 
We favor a preshock density somewhat higher than 
the canonical value of 0.2 cm""^. This may be related 
to the lag between the outer faint emission seen in Fig- 
ure 1 and the filament we observed. Thus while the 
soft rim is solid evidence for recent interaction with 
dense gas in the regions bright in X-rays and optical 
emission, its interpretation is still ambiguous in the 
Balmer-dominated filaments in the northern Cygnus 
Loop. 

4. Summary 

Far UV observations of the O VI lines in a non- 
radiative shock in the northern Cygnus Loop show 
that the oxygen and proton temperatures are fairly 
close to equilibrium; To/Tp < 2.5. This is con- 
sistent with the effective electron-ion equilibration 
in this shock inferred by Ghavamian et al. (2001) 
and effective ion-ion equilibration, as opposed to the 
strong preferential heating of oxygen measured in 
faster shocks in the solar wind. 

The O VI line profiles show clear evidence for 
resonance scattering, both in the profile shapes and 
in the intensity ratio of the A1032 and A1037 lines. 
While there is some imccrtainty as to the relative 
contributions of interstellar O VI opacity and opac- 
ity within the emitting sheet of gas, optical depths 
T1037 — 0.5 — 2 are required, with the highest optical 



depths in the bright filaments and at small velocities. 
The opacity effects reduce the total line intensity of 
the 1032 line by roughly a factor of 2. 

The line shapes and intensities make it possible to 
estimate the structure of the emitting sheet along the 
line of sight, and to infer a pre-shock density. These 
estimates can be checked against ROSAT data. We 
find a depth along the line of sight of 0.7 to 1.5 pc and 
a pre-shock density of 0.3-0.5 cm~^. A shock speed 
of 350 km s~^ at the upper end of the range allowed 
by the Ha profile (Ghavamian et al. 2001) is needed 
to match the FUSE and ROSAT data. The above pa- 
rameters support that conclusion to some extent, in 
that the pre-shock density is roughly twice the canon- 
ical value, but part of that difference results from the 
difference between the 440 pc distance to the Cygnus 
Loop and the 770 pc distance used in earlier analyses. 

Levenson et al. (1999) used the soft outer edge of 
the Cygnus Loop in ROSAT images to conclude that 
the blast wave is slowing down due to a very recent 
encounter with higher density material. The shock 
model shows that departure from ionization equilib- 
rium can contribute to the softness of the ROSAT 
spectrum in the outer 1-2' behind the Balmer line fil- 
aments, reducing the need for a density enchancement 
in this region. 
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Figure 1: Images of the northern Cygnus Loop 
Balmer line filaments in Ha (red) and [O III] A5007 
(green). North is at the top and East at the left. 
The image is 7.3' by 9.2'. An expanded view of the 4 
MDRS slit positions is shown in the inset. The 4 over- 
lapping LWRS exposures lie on the diffuse emission 
near the NW corner. Positions 1 through 4 progress 
from northeast to southwest. The 30" size of the 
LWRS aperture shows the scale. 



This 2-column preprint was prepared with the AAS lATJiiX 
macros v4.0. 
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Figure 2 - (a) O VI A1032 line profiles for the four 
MDRS slit positions, with the topmost profile corre- 
sponding to the leading filament. Profiles of Positions 
1, 2 and 3 are offset by 2.0, 1.5 and 1.0 to avoid over- 
lap, (b) Same for O VI A1037. 



Figure 3: O VI A1032 profile for the total nighttime 
exposure of all four positions for the LWRS. 
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Figure 4: HUT spectrum showing the wavelength 
range including C IV AA1548, 1550 and He II A1640. 
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Figure 5 -(a) The intensity ratio of Equation 2 for 
Position 1. Dotted line indicates velocity range where 
absorption by C II or H2 affects the ratio, (b) Same 
for Position 4. 
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Figure 6: Theoretical value of the intensity ratio R 
plotted against noa?. 
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Figure 7: Shape of the emitting sheet along the 
line of sight. Note that the X-scale is compressed 
by a factor of 5 compared with the Y scale. Dashed 
lines indicate the lines of sight for the 4 FUSE MDRS 
spectra. 
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Figure 8 -(a) Predicted O VI A1032 line profiles 
for the four MDRS spectra for a 350 km s"'^ shock 
with the shape of Figure 7 and a pre-shock density of 
0.5 cm-3. (8b) Same for O VI A1037. 
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Figure 9 -(a) Predicted intensity ratios correspond- 
ing R to Figure 5a for position 1. (b) Predicted in- 
tensity ratios corresponding to Figure 5b for position 
4. 
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Figure 10: O VI A1032 profile for the optically 

thin case and including the effects of interstellar and 
intrinsic resonance scattering. Interstellar reddening 
will reduce the fluxes by about a factor of 2.8. 




Figure 11: Predicted ROSAT PSPC count rates as 
a function of distance behind the Ha filament, along 
with measured values (+). 
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